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RAYSTC:  A COMPUTER  CODE  FOR  CALCULATING 
SINGLE  RAY-PATH  STATISTICS,  ASSUMING 
THE  GARRETT-MUNK  MODEL  OF  INTERNAL  WAVES 

I.  INTRODUCTION 

In  a recent  paper  [1]  Munk  and  Zachariasen  used  the  Rytov  approximation  to  calculate 
the  expressions  for  the  mean-square  acoustic  phase,  phase  rate,  log-intensity,  and  their  spectra 
for  a single  ray  in  an  ocean  possessing  internal-wave  induced  sound-speed  fluctuations.  They 
expressed  their  results  as  integrals  over  the  ray  path  of  quantities  associated  with  the  internal- 
wave  spectrum.  The  purpose  of  this  report  is  to  document  a computer  program  designed  to 
evaluate  these  integrals  for  the  case  in  which  the  internal-wave  spectrum  is  given  by  the 
Garrett-Munk  internal-wave  model  [2-4], 

In  Sec.  II  we  briefly  outline  the  theory  developed  by  Munk  and  Zachariasen  and  describe 
the  algorithm  used  in  the  computer  program.  Section  III  contains  a description  of  the  program. 
In  Appendix  A we  indicate  the  deck  assembly  and  list  in  a table  all  the  input  parameters.  In 
Appendix  B the  source  language  listing  is  given.  Appendix  C contains  the  results  of  a sample 
run  in  which  we  have  assumed  that  the  ray  path  is  quadratic  in  the  range  variable.  This  choice 
of  ray  path  was  dictated  solely  by  convenience,  and  is  not  necessarily  intended  to  be  representa- 
tive of  an  actual  physical  situation. 

II.  SOME  BASIC  RELATIONS 


According  to  the  general  Garrett  and  Munk  internal-wave  model  (2-4),  the  vertical  parti- 
cle displacement  { at  depth  z has  the  spectrum 


where 

F{(«,y;z)  - <{2(z)  >G(o*.  n(z))  H(j), 

(la) 

G(u,n)  » 0 for  «<«,  and  w > n, 

(lb) 

and 

f dw  G(u,n)  ■ 1, 

Jm, 

m 

(lc) 

£ H(J)  - 1. 

(Id) 

j-\ 


Here  w and  j are,  respectively,  the  internal-wave  radial  frequency  and  the  mode  number,  and 
<o,  is  the  inertial  frequency.  The  local  buoyancy  (Brunt-VaisalM)  frequency  has  the  form 

n (z)  - n0  exp(-z/£)  (2) 

in  an  exponentially  stratified  ocean.  The  mean-square  displacement  is  dependent  on  depth 
through  n(z),  i.e., 
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where  <{„>  is  the  mean-square  displacement  extrapolated  to  the  surface.  The  horizontal  and 
vertical  wave  numbers  are  given  by  the  dispersion  relations 

<4a) 


The  analysis  in  Ref.  1 is  based  on  a version  of  the  Garrett-Munk  model  [4]  GM75  1/2, 


where 


0 ; otherwise, 


The  mode  number  parameter  J*  is  typically  set  equal  to  3.  The  quantities  Sq,  Nh  are  dimen' 
sionless  normalization  constants  determined  from  Eqs.  (lc)  and  (Id): 


We  are  interested  in  the  pressure  received  at  a point  x - (Rm,„  0,  z ) in  the  ocean  due  to 
a source  located  at  x,  - (0,  0,  z,)  and  radiating  acoustic  energy  at  frequency  /.  Propagation  is, 
therefore,  along  the  x-axis  and  R m„  is  the  range.  (Since  we  consider  the  manner  in  which  the 
rms  values  accumulate  as  one  moves  along  the  horizontal  path  from  source  to  receiver,  we 
have  appended  the  subscript  max  to  indicate  R ^ is  the  maximum  range  of  interest.)  The 
acoustic  wavenumber  is  k„  - 2w f/c0  where,  typically,  c„  — 1.5  km-Hz. 

In  the  absence  of  internal  waves,  the  sound  speed  is  equal  to  the  (depth-dependent)  mean 
sound  speed: 

c(x,  /)  - c(z)  . (7) 

and  the  received  pressure  is 

Relp#(*)e-,"/>)  . (8) 
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In  Ref.  1 it  was  assumed  p0  could  be  approximated  using  ray  acoustics.  The  ray  path  of 
interest,  z(x)f,  satisfies  the  eikonal  equation 

-----  + K'(z(x))  - 0,  (9a) 


where 


and  the  end-point  conditions 


V'(z)  - - — — , 

2 dx  c(z) 


The  slope  of  the  path  is  given  by 


z (0)  — z, , z(RmJ-z. 


tan0(x)  — 


(The  coordinate  system  is  orientated  so  that  the  positive  z-direction  is  downward,  hence  the  ray 
angle  0(x)  is  positive  if  the  ray  is  directed  toward  the  bottom.) 

Also  of  interest  is  the  phase  curvature  function  [S]  A (x)  defined  by  the  expression 


M(x)]- 


where  f j 2 are  the  linearly  independent  solutions  of  the  differential  equation 

( d1  , r _ n 


+ V (z(x))|  f,,2  “ 0, 


that  satisfy  the  boundary  conditions 


r(2)  - -I  -Sl. 

W 2 dz 2 c(z)  ’ 


1 1(0)  - URmn)  - o, 


*|(*n»x)  - *j(0)  « 1. 


(Because  of  the  Wronskian  relation  obeyed  by  f | 2,  the  denominator  in  Eq.  (11)  is  actually 
independent  of  x.) 

When  internal  waves  are  present,  the  sound  speed  acquires  a small  fluctuation  Sc(x,  t ) 
which  is  related  to  the  internal-wave  vertical  displacement  by 


fWe  have  assumed  the  parabolic  approximation.  Consequently,  the  tec1*  term  which  appears  in  the  expressions  for 
the  mean-square  values  in  Ref.  1 has  been  replaced  by  unity  and  the  eikonal  equation  which  defines  the  ray  path  has  a 
slishtly  modified  form. 
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(constant)  r2(z)((x,  t) 


Therefore,  from  Eq.  (3), 
/ [8c(x,  i) 


is  the  mean-square  fractional  fluctuation  extrapolated  to  the  surface. 
With  the  presence  of  8c,  the  received  pressure  can  be  expressed  in  the  form 


where 


with  X\'2  real.  This  expression  can  be  rewritten  as 


U — <i>]  +2  7ri[<f>—  <^>) 


where  anticipating  the  use  of  the  Rytov  approximation  we  have  introduced  the  definitions 


(The  2 ir  is  included  here  because  we  have  chosen  to  measure  phase  and  phase-rate  statistics  in 
cycles  and  cycles-per-hour,  respectively.)  Mean-square  values  are  given  by  the  expressions 

<[*-  <*>12>  , (19a) 


Here  a dot  indicates  differentiation  with  respect  to  time  and  is  the  rms  intensity  in  decibels. 
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Using  the  Rytov  approximation,  Munk  and  Zachariasen  [1]  calculated  the  contribution  to 
these  mean-square  values  from  a single  ray  z(x).  They  found f 

4>L  “ J*0Rm“  * (z(x))J  L,  («(x).  z(x))  , (20a) 

(zU))|  Vp(0(x),  z(x))  , (20b) 

<IM>- !<xi)  <2oci 


In  these  expressions  we  have  introduced 


Lp(9.  z)  - L(0)fx  tan0(z)  , 


uo)  = <r‘>  I— U 

w 


1 , A2  . (A2+l)1/2  + l 

A2  + 1 2(A2  + 1)3/2  (A2+l)'/2-l  ’ 


V,C.z)  - ± <r’>  *,n0B  In  tan  «J. 


JL  f(A)  = JL  Ai  _ I 1 . (A2-H)i/2-H 

/2<  “ 2/n  4 2 (A2+l)*'2  (A2-!-!)1'2-! 


L\(.z ) ■ 


1 [<r'>  (»;.-i) 


a 

«(r)  B 


The  quantity  <j~'>  is  the  average  of  the  reciprocal  of  the  mode  number  over  the  internal- 
wave  spectrum*: 

<r,>-"”|77T7T'  (14> 


fSec  footnote  on  p.  J. 

’Properly,  one  should  write  /j  as  t function  of  both  4 and  n/u,. 

*Uaing  a technique  which  we  believe  to  be  accurate  to  within  I part  in  10*.  we  And  <J~X  > - 0.7308,  0.6240,  0 4890 
0.4001,  0.3404,  0.2971  for  J.  - 0,  1,  2,  3,  4,  5,  respectively. 
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All  the  other  quantities  in  Eqs.  (20)-(23)  have  been  previously  defined.  The  function  Lp(9,z)  is 
the  sound-speed  correlation  length  measured  at  depth  z along  a line  inclined  at  an  angle  0 with 
the  horizontal,  and  Lv(z)  is  the  vertical  correlation  length  at  depth  z.  The  functions  /i(A)  and 
/j(A)  are  defined  in  such  a way  as  to  be  equal  to  unity  when  A is  equal  to  zero.  The  program 
assumes  n„  and  w,  are  given  in  cycles-per-hour.  With  this  choice,  4>tm,  is  measured  in  cycles 
and  is  measured  in  cycles-per-hour.  In  the  program  the  integrals  Eqs.  (20)  are  calculated 
using  the  trapezoidal  rule.  The  range  of  propagation  is  broken  up  into  NR  slabs  of  equal  width, 
and  the  integrands  in  Eqs.  (20)  are  assumed  to  be  linear  within  each  slab.  For  all  the  cases  we 
considered,  this  simple  algorithm  gave  completely  satisfactory  results. 

III.  PROGRAM  DESCRIPTION 

The  program  was  written  for  use  on  the  Texas  Instruments  Advanced  Scientific  Computer 
in  the  machine-specific  language  TI-ASC  FORTRAN.  Since  it  is  a complete  program  rather 
than  a subroutine,  its  use  does  not  require  a calling  sequence. 

The  total  length  of  the  program  is  0001A600.  The  system  reserves  an  additional  8K 
words  of  central  memory  for  I/O  buffers,  etc.  It  requires  no  temporary  storage  and  does  not 
use  common  blocks.  The  program  will  compile  on  the  NX  compiler  at  the  K level  in  1.60  s. 
The  execution  time  varies.  For  the  sample  computer  run  recorded  in  Appendix  C,  the  total 
central  processor  time  was  5.11  s and  the  plotter  time  was  14  min. 

The  program  uses  the  following  external  routines:  ABS,  ALOG,  ATAN,  EXP,  FLOAT. 
INT,  SQRT,  TAN,  ORIGIN,  NXAXIS,  NYAXIS,  LETTER,  NUMBER,  PLOTS,  PLOT, 
ENDPLT,  RJTOP. 

The  required  input  data  are  listed  in  Table  A1  of  Appendix  A.  They  naturally  fall  within 
four  categories:  acoustic  parameters,  internal-wave  parameters,  output-option  parameters,  and 
ray  characteristics.  The  acoustic,  internal-wave,  and  output-option  parameters  are  entered  on 
three  separate  cards  which  together  compose  an  input  file  embedded  in  the  job  input  stream  by 
means  of  a START/STOP  statement  pair.  This  file  has  the  standard  Fortran  access  name 
FT05F001  and  hence  is  read  on  logical  unit  number  5.  Distances  are  given  in  meters,  the 
acoustic  frequency  is  given  in  Hertz,  and  the  frequencies  associated  with  the  internal-wave 
model  are  given  in  cycles-per-hour. 

The  present  version  of  the  program  does  not  calculate  the  ray  path  nor  the  phase  curva- 
ture function.  The  depth  of  the  ray  path  as  a function  of  horizontal  path  length  is  input  by 
means  of  a card  file  specified  by  a START/STOP  pair  and  having  the  access  name  FT08F001 
(logical  unit  number  8).  This  file  consists  of  the  depth  of  the  ray  path  at  NR  + 1 equally-spaced 
range  points  extending  from  the  source  position  to  Rmtx  (RMAX).  These  values  are  stored  in 
the  array  RAY  in  such  a way  that  RAY(I)  (1*1,  ....  NR+1)  is  the  depth  of  the  ray  path  (in 
meters)  at  a horizontal  distance  (1-1)  RMAX/NR  from  the  source.  The  maximum  value  of  NR 
allowed  by  the  program  is  7000. 

It  is  not  the  phase  curvature  function  which  is  required  as  input  but  rather  the  absolute 
value  of  its  reciprocal.  Just  as  with  the  ray  path,  NR+1  values  of  this  function  (in  meters) 
compose  a card  file  specified  by  a START/STOP  pair.  The  file  has  the  access  name  FT09F001 
(logical  unit  9).  These  values  are  stored  in  the  array  ABRECA  so  that  the  absolute  value  at  a 
horizontal  distance  (I- 1)  RMAX/NR  is  ABRECA (I)  where  1*1, ....  NR  + 1. 
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On  encountering  an  error  condition  in  the  transfer  of  input  data,  the  program  will  write 
out  the  status  code  of  the  error  message  and  terminate  the  job. 


The  output  from  the  program  can  logically  be  divided  into  three  categories;  /.  input 
parameters,  derived  quantities  of  secondary  importance,  and  results,  /'/.  optional  tables  of  vari- 
ous arrays,  and  iii.  optional  plots.  In  the  following  three  paragraphs  we  will  describe  these 
categories.  All  output  is  written  on  the  standard  printer  (logical  unit  6).  The  standard  access 
name  FT59F001  (logical  unit  S9)  is  assigned  to  the  plotter  output. 


The  program  lists  the  input  parameters  from  the  first  two  cards  of  the  input  file  FT05F001 
(Table  Al)  together  with  their  Fortran  names  and  units.  In  addition,  the  source  and  receiver 
depths  are  listed.  These  depths  are  the  values  of  the  first  and  last  elements  of  the  ray-path 
array  RAY.  The  program  lists  the  values  calculated  for  the  range  increment  A r — 


X/NR  (DELR),  the  wave  number  k0  (WV),  the  value  for  1(0)  (LZERO),  calculated  from 


Eq.  (21b),  and  the  value  for  <j~x > (AVE).  The  program  then  lists  the  values  obtained  for 
the  root-mean-square  phase,  phase-rate,  and  intensity  fluctuations  together  with  the  errors  asso- 
ciated with  the  use  of  the  trapezoidal  rule. 


The  user  has  the  option  of  listing  in  tables  the  values  of  the  elements  of  six  arrays  for 
selected  values  of  the  indices.  These  arrays  are: 


1.  RANGE  — gives  values  in  meters  for  the  horizontal  path  length  at  NR  + 1 evenly  spaced 
points  along  the  horizontal  path  of  propagation. 


2.  RAY  — contains  values  for  the  depth  of  the  ray  path. 


3.  ANGLE  — contains  the  NR  values  for  the  ray’s  grazing  angle,  in  degrees,  calculated 
using  a finite  difference  approximation  to  the  derivative  of  the  ray  path. 


4.  ABRECA  — contains  values  for  the  magnitude  of  the  phase  curvature  function.  This 
array  is  labelled  |(1/A)|  in  the  Appendix  C table. 


5.  PHI  — contains  values  for  the  calculated  rms  phase  fluctuation  along  the  ray  path  as  a 
function  of  horizontal  distance  from  the  source. 


PHIDOT  — contains  values  of  the  rms  phase-rate  fluctuation  as  a function  of  distance 
from  the  source. 


The  extent  to  which  the  tables  are  constructed  is  determined  by  the  output-option  parameters 
IARR  and  NPRNT  (see  Table  Al). 


Depending  on  the  values  assigned  the  output-option  parameters  NPLT  and  IPLT(I),  I — 

1 7,  the  program  will  construct  up  to  seven  plots.  Examples  of  these  plots  are  given  in 

Appendix  C and  Table  Al  contains  brief  descriptions. 
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Appendix  A 

DECK  ASSEMBLY 

Card 

Number 

Description 

JSL  Statement  Format* 

1 

JOB 

/bJOBbName,  Acct.  no..  User  code,  OPT  — (C,  R,  D,  T) 

2 

LIMIT 

/bLIMITbMIN  “ 1,  BAND  - 25 

3 

PLOT  FILE 
DESCRIPTION 

/bFDbFT59F001,  FORG-PS,  RCFM-U.BKSZ-4000, BAND- 1/10/1 

4 

FORTRAN 

/bFTNblN  - SDECK,  FTVERS  - NX,  FTNOPT  - (K,U) 

5 

LINK 

/bLNK 

6 

EXECUTE 

/bFXQTbOPT  - (I, A) 

7 

PLOT  OUTPUT 

/bFOSYSbFT 59F001 , TYPE  - PLOT 

8 

START  SOURCE 

/bSTARTbACNM  - SDECK 

DECK 

Fortran  source  deck 

9 

STOP  SOURCE 
DECK 

/bSTOP 

10 

START  INPUT 

/bSTARTbACNM  - FT05F001 

PARAMETER 

Input  parameter  cards 

11 

STOP  INPUT 
PARAMETER 

/bSTOP 

12 

START  RAY 

/bSTARTbACNM  - FT08F001 

ray  path  data  cards 

13 

STOP  RAY 

/bSTOP 

14 

START  PHASE 

/bSTARTbACNM  - FT09F001 

CURVATURE 

phase  curvature  data  cards 

15 

STOP  PHASE 
CURVATURE 

/bSTOP 

16 

END  OF  JOB 

/bEOJ 

•The  lower  case  letter  “b"  represents  a (necessary)  blank. 
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Table  A1  — Input  Data  For  RAYSTC 

FILE  ACCESS  NAME  FTOSFOOl 

Format  (F7.2,  F8.3,  F8.0,  15,  F5.0) 

CZERO 

Reference  sound  speed  c,  (m/sec) 

FREQ 

Acoustic  frequency  / (Hz) 

RMAX 

Maximum  range  R (m) 

NR 

Number  of  range  steps 

DEPTH 

Ocean  depth  D (m) 

CARD  2 

Format  (F5.0,  F5.2,  E9.2,  F6.4,  15) 

B Stratification  scale  B (m) 

BVZERO  Buoyancy  frequency,  extrapolated  to  the  surface,  fl„(cph) 

DELC  Root-mean-square  value  of  the  tract 'ona^sound  speed  fluctuation, 
extrapolated  to  the  surface,  <(fir/r)^> 

FREQIN  Inertial  frequency  ».  (cph) 

JSTAR  Mode  number  parameter  j. 

CARD  3 Format  (11,  14,  811) 

IARR  Table  output  flag 

— 0;  do  not  print  tables 

— 1;  print  tables  of  the  arrays  RANGE,  RAY,  ANGLE,  ABRECA,  PHI,  and  PHIDOT 

NPRNT  Table  output  count;  the  tables  will  be  composed 

of  the  first,  last  and  every  NPRNTth  element  of  the  arrays. 

NPLT  Plot  output  flag 

— 0;  construct  no  plots 

— 1;  construct  one  or  more  of  the  7 possible  plots 

IPLT(I)  Individual  plot  flag 

— 0;  do  not  construct  plot 

— 1;  construct  plot 

IPLT(l)  — Plot  of  L,(9,z)  vs  9 for  various  values  of 
z (Fig.  1 of  the  sample  output  of  App.  C) 

IPLT(2)  — plot  of  the  ray  angle  9 vs  horizontal  path 
length  (Fig.  2 of  App.  C) 

IPLT(3)  — plot  of  the  ray  path  (Fig.  3 of  App.  C) 


IPLT(4)  — plot  of  the  absolute  value  of  the  reciprocal  of  the  phase 
curvature  function  (Fig.  4 of  App.  C) 

IPLT(S)  — plot  of  the  rms  phase  fluctuation  vs  horizontal  path 
length  (Fig.  S of  App.  C) 

IPLT(6)  — plot  of  the  rms  phase-rate  fluctuation  vs  horizontal 
path  length  <Fig.  6 of  App.  C). 

IPLT(7)  — plot  of  t^(#,r)  vs  9 for  various  values  of  z (Fig.  7 of  App.  C) 

PILE  ACCESS  NAME  FTMFM1 

This  file  contains  a number  of  cards,  format  (10F8.1),  listing  the  depth  in  meters  of  the  ray  path  at  NR+1  evenly 
spaced  points  along  the  horizontal  direction  of  propagation. 

FILE  ACCESS  NAME  FTHFM1 

This  card  file  is  analogous  to  FT08F00I  but  contains  values  in  meters  of  the  reciprocal  of  the  phase  curvature  function. 


10 


Appendix  B 

SOURCE  LANGUAGE  LISTING 


FR06RAN  RATSTC 

DIMENSION  RAT(TOtl). ABAECA(TQOl). ATE J(4).  SLOPE (7 000 )»ANSIE (7000). 
1PRISQ<2. 7001 >.DOTSQ(2. 70)1 >.I0TASQ(2. 7001 >.FHIC7001>*PNIDOTC7001>» 
2PLTAROOOO).  IFLK7) 

DIMENSION  XNAX<7>.TNIN(7>*THAX(7>tTINC<7).NDIST<7).BCXl(4>*0CX2<4) 
1.0CT1(4).BCT2(4).BCT3C2).BCT4(3).BCT9(3).BCT6(3).BCT7(4) 

DATA  AVEJ/0. 7300. 0.4240. 0.4090. 0.4001 *0.3404. 0.2970/ 

DATA  DCI/2NKN/ 

OATA  HOIST/-!. 1. -1.3. 3.3. ~1/ 

OATA  BCX1/4HAN6L.4HE(DE»4H6REE.4HS)  / 

OATA  SCX2/4MR4NS.4HE(KN.4H>  / 

DATA  RCT1/4NL<P>.4H<N-C*4HTCLE.4H442)/ 

OATA  RCT2/4NANSL.4NE(0E*4HSREE.4HS>  / 

OATA  DCTD/4NAATC*4NN>  ‘ / 

DATA  DC 74/4HA0RE . 4NCA( X * ANN > / 

DATA  0CT9/4NPNIC.4NCTCL.4NES)  / 

DATA  DCT4/4NFNID.4N0T(C.4HHP)  / 

DATA  0CTT/4NV<P).4N<N-C*4NFN4».4N2>  • 

CALL  MiT# 

INTESEA  ONE 

REALM  LIEN0.L9.LVIEA0* IOTASO*IOTOT.NU 

IERO-0.0 

WRITE  (4.  97) 

MRNATCINl) 


RE A0<9«1.ERR«4.STV*ISTAT)CZER0.FREQ.RNAX.NR.DEFTN 
FOANAT(F7.2.FS.3.FO.O. I9.F9.0) 

REAO(9.2.ERAM.STV«ISTAT)B.DVZERO.DELC.FREQIN.  JSTAR 
FORM AT (F9.0. FS.2.E 9.2*  F4 . 4. 19  > 

READ(9.3.EAA«4.STV»ISTAT)IAAA.NFRNT.NFLT.CIFLT(I)*I*1*T) 

FORNATdl.I4.OII) 

SO  TO  A 

WRITE <4.7)11 TAT 
F0RNATC1N  .SNERROR.IlO) 

STOF 
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••39 

••«• 

••41 

••42 

••43 

••44 

••49 

••44 


•ETAa»VZERO/FREQIN 
AVE-AVEJEJSTAR  * t> 

UiH><W(»KTM)/(nM4) 

UMM(((n*D(lC)  )M|)«UflM|U 

BCLT«-«(UV*OCLC>  )M2)*VZER4*DELA 
LVZERNCB/RISTARlFSQRWAVERCRISTAR-m/RI) 
«»>(((ru*U)/tiK(ti.))Mt)  /(w*LncRMi«ER«) 


I 


H 


(I 

* i 


C REM  IN  RRT  RATH  AHO  ABSOLUTE  VALVE  OF  RECIFROCAL  OF  FHASE  CURVATURE  FUNCTION 

M4T  REAOOtt  l*ERRa4f  STV«ISTAT>(RAT<I>»IaI*NRFl) 

••41  REAOO.Ll.ERRX.STV-ISTATHAaRECAID.I-l.NRRI) 

••49  U F0RHATO •••«!) 

MS*  ISCR-RAT(l) 

••SI'  ZREC-RATCNRFt) 

C CALCULATE  SLOFE  AHO  LOCAL  ANCLE  RAT  RATH  HARES  HITN  RESRECT  T9  THE  HORIZONTAL 

••92  ••  S«  I~t.HR 

••S3  SL0RE<I>-CRATCI*1)-RAVCI»/0ELR 

••94  AN6LE<I3aATAN(SL4RECI3) 

••99  9*  CONTINUE 


F 

1 1 


| 


i 

[. 


C CALCULATE  RN1S9  OOTSQ  AND  IOTASQ  ARRATS 


••9* 

N9T 

MS* 

••99 

MM 

•••1 

••42 

••43 

••44 

••49 

••44 


N It  J*I#I 
RHISQCJi  !)>••• 
OOTSRCJt !>■••• 
IOTASQCJ »!)-•.• 

••  l(  I- I .HR 

T1 *C ANNAOR A T«*4- 1 ) 

TZ-EIRCSatFTI) 

T3»E*F<Tl) 

T9a(IR(2.MTl) 

X4aRETA*T3 

IaI4oASS(SL0RECI3) 


C CALCULATE  El  AHO  E2  FUNCTIONS  NAVINC  CALCULATED  X 


NAT 


N49 

NT* 

NT1 

NT2 

NTS 


KSNXOX 

IF(XSQ.NE.ZER0)60  TO  it 

n»i.t 

E2aI*« 

64  TO  32 
34  ItMSMlJ 
X2aS0RT(II) 


12 


■mi 


NRL  REPORT  8304 


X3>X2**3 

X3-XSQ/X3 

FXIM.I/I1 

FI2<«IW(X2«1.I) 

FX3-XS*FX2 

FX4-.9*XS««l9t(X$Q) 

FX»FXl  4F  X3-FX4 
FXMUK(X4) 

FXt>.i«tlK  (XSQ/4) 

FXT«(U.0C(X2»t.l)/X2)-(AL0C(XSQ)/(2*X2)) 

F2>(FX5-FX4-FXT)/FXS 


S-M.FIM«T2*FX 

U-NU*T9»N0«ECA<I»J-1) 

PNISQC  J*  IMWHISQC  J»  «♦$ 

I0T«SQ(J»I»1  >-IOTASQ(J»D»CS*U> 

DfTSQC  J*  1*1  )«OOTIQ(  J*  I )«COELTA«T2*AL06(X4)*F2) 
CONTINUE 


Of  33  I-l.NNM 

Dl*0.S«<FNISQ<lt»»PNISQ<2tI» 

02-3  .StCDOTSOCl  1 1 )*O0TSQ(2.  IM 
FNI(I>sSQRTC01> 

PH  ID  OT  C I >“S9t T <02 ) 

CONTINUE 

FNIT0T*FNI(NRF1 > 

OOTTOT«FNIOfT(NRFl) 

IOTfT-SONT(f.S«(IOT«SQCi.NRFI>»  I9T«SQ(2*NRF1)» 
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0112 

92 

0113 

300 

011A 

otto 

301 

011A 

OUT 

302 

0110 

0110 

303 

0120 

0121 

30  A 

0122 

0129 

309 

012A 

0120 

30A 

0124 

012T 

307 

0120 

0120 

300 

0130 

0131 

309 

0132 

0139 

310 

0134 

0130 

311 

0134 

0137 

312 

0130 

0130 

313 

01  AO 
01A1 

31A 

01A2 

01A9 

319 

01AA 

01AO 

31A 

01AA 

01AT 

317 

01A0 

OlAO 

310 

0100 

0101 

319 

0192 

0193 

320 

MS  SUM  INC  THE  (1M(TT<MIK  MODEL 
HRITECAt  301) 

FORNATClHOtZSHINPUT  ACOUSTIC  PAR 

HRITECAt  302)C2ER0 

FORNATC1H  t33HREFERENCE  SOU HO  V 

HRITECAt  30 3 ) FEE Q 

FORNATCtN  t33N ACOUSTIC  FREQUCNC1 

HRITECA.30A)RNAR 

FORM  At  OH  • 33HNAX  RANCE  (ROAR) 

HRITECAt  309)NR 

FORNATCtN  t)3HNUH0ER  OF  RANCE  SI 

HRITECAt  30A)0EPTH 

FOR NATC1N  t33NOCEAN  OEFTN  COEFTH 

HRITECAt 307)ZSCR 

FORNATC1N  .33HSOURCE  OEFTN  CISC! 

HRITECAt  30DZREC 

FOR  RATON  tSSNRECEIVER  DEPTH  CZR 
HRITECAt 

FORRATONOt 30NINPUT  INTERNAL  HAT 
HRITECAt 910)0 


HRITECAt  91D0VZER0 
FORNATCtN  tAZNERTRA 
IPN) 

HRITECAt  91Z)0ELC 
FORNATCtN  tAtNERTRA 
HRITECAt 91 9IFREQIN 
FORNATCtN  tAZNINERT 
1FH) 

HRITECAt 31 A) JSTAR 
FORNATCtN  tAZNJSTAA 
HRITECAt 919) 
FORRATClNOtZlNCALCU 
HRITECAt 91 A)0CLR 
FORNATC1N0 tZIHRANCE 
HRITECAt 31T3HV 
FORNATCtN  tZlNHAVEM 
HRITECA.SIOILZERO 
FORNATCtN  tZINLZERfl 
HRITECAt 31«) ATE 
FORNATCtN  tZtNAVE  < 
HRITECAt  320) 


INTERNAL  HAVES/) 


•PT.ZtAN  N/SEC) 
tPO.StAN  HERTZ) 
tFO.OtZN  N) 
tIS) 

tFS.OtZN  N> 
tFA.ltZN  N) 
tFA.ltZN  N) 

- tFO.OtZN  N 

RO)  - tFS.ZtAH  C 

DELC)  - tEO.Z) 

» tFA.AtAN  C 

- til) 


40 ER  CHV)  ■ ,FA.At3N  /N) 

* tFO.St  3H  N) 
t/J  CAVE)  • tFA.A) 
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WRITECA.321 IPHITOT 

F0RNATC1NR.4 1HRNS  PHASE  FLUCTUATION  CPHITOT) 
1VCLE) 

WRITECA. 322 JO0TT0T 

FORNATCIN  .41HRNS  PHASE  RATE  FLUCTUATION  COOTTOT) 
IPN) 

WRITE C6. 323 )IOTOT 

FORNATCIN. .41NRNS  INTENSITY  FLUCTUATION  (IOTOT) 

1 00) 

WRITECA. 324 )EPHI 

FORNATC1  HA.SAHERROR  FOR  PHASE  (EPHI)  - »F 
14.2. AN  PERCENT! 

WRITE <«• 323 HOOT 

FORNATCIN  tSRHERROR  FOR  PHASE  RATE  CEOOT)  • .1 

14.2. OH  PERCENT) 

WRITE  CA.32A)EI0T 

FORNATCIN  .SOHERROR  FOil  INTENSITY  CEIOT)  - i 
IF*. 2, AH  PERCENT) 


WRITE CA.32T) 

F0RNATC1 HO .AHT ARLES) 

WRITECA.32A) 

FORHAT  C I HA  . ASHINOE  X RANCE-H  RAYCD-N  ANGLE  (t-l)-RAO  1C 
U/A)ICI)-H  PNICD-CTCLE  PHIOOTCD-CPN) 

ONE- I. A 

WRITECAt329)0NEf  ZERO.RATCD*  A0RECAC1)  .PMICD.PHIOOTCl) 
F0RRATClNO*ISt4X.FA«A*3X.FA.lt24X.Fi.I.TX.FA.S.4X.F9.S) 

SNR-CK-D-DELR 

NRITECA.33A)K»SNR»RAYCK)»ANCLECK-l)t AORECACK)  ,PHICR),FHIOOTCA) 
FORNATCIN  .IS.4X.F0.A.3X.F0.1 .4X.FA.S.I2X.FA.1.TX.F9.9, AX.FA.S) 
CONTINUE 

IPCR.EO.NRPDG0  TO  AO 
SNR-NRPOELR 

WRITECA. 33A)NRP1. SNR. RATCNRPI). ANGLE CHAP 1~1).A0RECACNRP1).PHICNRP1 
D.PHIOOTCNRP1) 


M.  L.  BLODGETT 


l!M  IF  IK  ST- 1 


C CALCULATE  A INI  AIM  AND  AAXIAUN  VALUES  ATTAINED  AT  FLATS 


DIES 

DID* 

DIET 

EISA 

• IE9 
fl«S 
Ei*l 
1192 
E193 
E194 
E19S 

• 19* 
E19T 
E19E 

• 199 
E2EE 
E2E1 
•292 
•293 
•2E4 
•299 
E2E* 
•2ET 
•2EE 
•2E9 
•2  IE 


BE  999  K-ltT 
TNIN<R3-B.B 

IF<K<EE<23  TAINCK3— 2».» 

IF<K<EE<I3  TAINCK 3-AEFTN 
999  CENTINUE 

VAAX<13-1E.«»INTK1.1*LZEAE3/1I.»3  ♦ 1E.E 
TNAK23-2A.* 

VAAX<33>t.» 

TEF-ABRECAU3 
•A  99T  K-ltNA 
CEAF - AAR  EC A<  R ♦ 1 3 
IF<CEAF<CT<TEF3  TAF-CENF 
99T  CENTINUE 

TAAXC43-1.1KTEF/1MS.3 
VAAXC93-1<1*FNITET 
TASK*)- 1.1*  OAT  TAT 

VAAXCT>-1E.»*INT<VZEAE*<ALECCEETA>*1>/10.0>  ♦ 1E.E 
DA  99*  K-ltT 

TINC<K3-<THAX<K3-THIN<K33/1D.> 

IF<K<NE<23  CE  TE  99* 

TINCC23-9.E 
99*  CENTINUE 

•A  99E  K-ltT 
KNAKCK3-AAAK/1EEI.E 
tF«K.EE.13.EA.<K.EE.T33KNAX<K3-lE.E 
99D  CENTINUE 


C MAN  FLETS 


•211  EE  TEA  K-ltT 

•212  IF<IFLT<R3<EE<A3  CE  TE  TM 

C BEAN  X ANB  T AXES 


•219 
•214 
•219 
•21* 
•2nT 
•21 A 
•219 
•22A 
•221 


•AAX>XAAX<K9 

ZAIN-TNIMKT 

ZAAX-TAAKK3 

XIA-ENAX/1E<E 

TIN-TINC<K9 

NX— I 

NT-NBICTCK3 

VSCALC-C2AAI-ZNIA3/1*.* 
X SCALE-KIN 
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miM 

ff(VIUT.«.l)  60  T8  19 
IPIRST-8 

CALL  00I6INC2.0  »•> 

I CALL  nUU(liMiinilUlill.i-.ITt.8Tin) 
V((.n.l.«4.N.I)  60  T0  13 
CALL  LtTTtAC««t.-.9..29*0CX2.0.12) 

I CALL  NVAXISC0.9ttNIN.UN.ZNAX.  1M  t-.IT. .0T.NV) 
CALL  MMKK'.lHt-.lllt.lMintH.taT) 

CALL  NXAXIS<l«19»0.9.XIN«9NAI.10*9ia07.9.NX> 

CALL  NVAXZSC10.0.9.2NZN. VINiZNAX.10«0i*67.0.NY) 

60  T0C001.002.003.604.609.006.0073.K 

P10T  1 LCP3VNS  0 


ILtL-1 

CALL  LfVVtRC3.9.~.99..29.BCX1.0.M> 
CALL  LEVVfAC-.979.3.t.*29.0CV1.90..14> 
•A0-PI/100. 

X AM -0.Z/X  SCALE 
00  A3  >IiHHiI9I 
IPCJ.EQ. 1299.0A.J.EO.1T9I)  60  T0  63 
IL6L-IL6L  ♦ 9 

ipiuo 

XAXIS«-XAM 

n 

00  <9  1-1.91 
TNE7A-(I-1>*.2*AA0 
X-7ANCTNETA)  Ml 


J 
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•249 

UHUIIMH) 

•244 

tFCLF.LT.9.99LF»9.9 

• 249 

VAXIS-LP/VSCALC 

•244 

CALL  FLETCXAXIStt  AXIS*  HEN) 

•24T 

IPEN-2 

•244 

IFCILOL.EQ. I)  II  T9  99 

•249 

U T9  29 

•2T9 

99 

2MN-J/1999. 

•2T1 

TAX. TAXIS  ♦ .93 

•2T2 

CALL  NUNIEICMIIS  ♦ .99.TAX  ..OT.ZNUNtE.Z) 

•2T9 

CALL  LETTEMIAIIS  ♦ .33. TAX  ,,.9T»9CI.9i2> 

•2T4 

CALL  FLETCXAXIStTAXISt 9) 

•2  T9 

29 

CMTINUE 

•2T4 

49 

CMTINUE 

•2TT 

M T9  41 

C FL9T  2 ANCLE (A)  VIS  1 

• 2T9 

••2 

XAXIS>9 

•2T9 

CALL  LETTEM-.STSt 3.9. .29* 0CT2.99.tl4) 

•2M 

09  21  HltNI 

•291 

XAXIS-CXAXIS*XAOO> 

•292 

0C6-109./PI 

•21) 

0C6IEE-AMLECI)9BES 

•294 

IFC0CCIEE.ST.29.996a  T9  12 

•299 

1FC0EMCC.lt  .-29.990E6IEE— 29.9 

•294 

60  T9  23 

929T 

12 

0C6IEC>29.9 

9299 

2) 

TAIIS-C0EMEE-ZNIN9/TSCALE 

9299 

CALL  PLETCXAXIStTAXIStIPCN) 

•299 

IFEN>2 

•291 

21 

CMTINUE 

•292 

60  T9  41 

C HIT  9 MTU)  VIS  I 


1299  119  MIIS— 1*00 

•294  CALL  LETTEIC-.9T9t4.lt. 29.IC79.99. .9) 

•299  II  2«  I-ltNIFl 

•294  IAXIS-CXAIIS«IA00> 

•29T  IFCAATCI>.LT.9.9>AATCI>>9.9 

•299  IPCIAVCI9.6T.OCFTH9IAVCI9-OCFTN 

•299  VAXIS.CAAVCIWNIN)/T$CALE 

9199  CALL  FL9T(  XAXIS  tT  AXIS  t HEN) 

• Ml  IPBN-2 

•9*2  29  CENTIME 

•999  II  T«  41 
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C HIT  4 UKCUI)  VRS  R 


•9*4 
• MS 
•Ml 
•MT 
•Mt 
•MS 
•111 
Oil 
•SIX 


••4 


4»S 


CALL  LETTEAC-.379t4.l(.29tSCT4(M.tl2) 
II  MS  I>1.NRP1 
XAXIS-XAXIS  ♦ MM 

TAXIS-CCAMECACD/IMI.WRIO/TSCALE 
CALL  PL0TCXASIS*7AXIStIPER) 

HIM 
C0NTINM 
M 70  41 


C PLOT  S MICA)  MS  I 


•SIS 

•S14 

•SIS 

•S14 

•SIT 

• SIS 
•SIS 

• SM 
•SX1 
•SX2 


MS 


22 


XAX1S—XAM 

CALL  LETTERC-.3?S*4.«*.2StlCT9tM.*12> 
••  22  I«1»HRP1 
X AXI  SS(X  AXIS  AXMS  > 
IPCPNICD.LT.OOPNICDM.I 
VAXIS-CPNKIWNIIO/TSCAIE 
CALL  PL0T<XAXIS«7AXIS*IPEH> 

1PER-2 
CMTIME 
40  T0  41 


C PLOT  4 PNIS0TCR)  MS  R 


•923 

•324 

• 92S 

• 924 
027 
•3  2* 

• 324 

• SM 
•SSI 

• 332 


Ml 


49 


XAXIS— XAM 

CALL  LETTER<-.3TS.4.0.2StBCT4«M.*12> 

00  IS  IflfRRPl 

KAIIS-CXAXIS«XAM> 

IPCPN100T<n.LT.I.OPNI00TC»«».S 

7AXISs(PNIO0TCI)-2AIN)/?SCALE 

CALL  PL0T<XAXIS( TAXIS* IPER) 

IPE0-2 

C0HTIMHE 


C PLOT  7 TCP)  MS  0 


• 933 

• 994 

• S9S 
034 

• 997 
•9M 


M7 


•A0-9.141SS24S/1M. 

CALL  LETTER<3.9a>.S«i.2St>Cll(«>14> 
CALL  LITTER<-.379*3.»*.2StK77tM.»ll) 
XA00M.2/XKALE 
M O JMt*M«2S« 
IPCJ*E«.12S».M.J.Ea.l7S»>  40  70  «3 
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• 33* 

IFEN-3 

• 341 

XAX1S— XAOO 

• 341 

T|»CXF(6ARIM*J)»B£TA 

•342 

T2-ALOS  (71) 

• 343 

••  •*  1*1*91 

• 344 

THETA*CI-l)*'.2*«AO 

• 349 

X*TAN(THETA>*T1 

C CALCULATE  F2  FUNCTION  HAVIN6  CALCULA7E0  X 

• 344 

xso*x*x 

• 347 

IF(XSQ.NE.IEft«>SO  70  9» 

• 340 

F2-1.0 

• 349 

SO  TO  91 

• 394 

9t 

Xl-XSQ  ♦ !•• 

•391 

K2-SMTCX1) 

•392 

FXS-AL0SCT1) 

•393 

FX4-.9*AL06  (XSQ/4) 

• 394 

FXT-<AL«S<X2*t.l)/X2)-<ALOC(XSQ>/C2*X2» 

•399 

F2*(FX9-FX4-FX7>/FXS 

• 394 

91 

XAXIS-XAII S ♦ XAOO 

• 397 

VFaVZER«*T2*F2 

• 390 

IFCVF.LT.0.»3VF*».» 

• 399 

TAXIS-VF/TSCALE 

•349 

• 341 

IFEN-2 

• 342 

IF (t. ME. 4)  SO  TO  OS 

• 343 

nWNJ/lNI. 

• 344 

TAXaTAXIS  ♦ .14 

• 349 

CALL  NUNMAfXAXIS  ♦ .•9*TAX*.»T.ZNUN.».2> 

•344 

CALL  LCTTER(XAXIS  • .33*TAX,.»T.0CI,t,2) 

• 347 

CALL  FL0T(XAXIS.TAXIS*3> 

• 34* 

•9 

CONTINUE 

•349 

•3 

CONTINUE 

•3T« 
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CALL  MISW12.»*»> 

•371 

7M 

CMT1NUE 

• 372 

CALL  ENOFLT 

• 373 

91 

STOF 

• 374 

ENO 

Appendix  C 
SAMPLE  RUN 


BLOOM 1 •OPT-(CtRtOtT) 


JOB  MARILYN  BLODGETT • 

LIMIT  HIN  = 1,BAND*<!5 
FO  FT59FCU1»F1RG*P  S»RCFM=U  »6KSZ **005 (BAND  = 1/10/1 
FTN  IN=SOECK«FTVERS=NX»FTNOPT=CKtU) 

LNK 

FXQT  <JPT*(I,A) 

FOSYS  FT59FO0 1# T YPE=PLOT 
START  ACNN-SOECK 


FORTRAN  SOURCE  DECK 


/ STOP 
/ START 
1500.00 
1000.  3. 

1 1011! 

/ STOP 

/ START  ACNN*FT08F001 


ACNM-FT05F9C1 

11/0.000  10000. 

CO  Q»49F-i 3C .)4  20 
11111 


1C04040. 

3 


1000. c 

960.4 

921.6 

003.6 

046.4 

010.9 

774.4 

739.6 

705.6 

672.4 

040.0 

600.4 

577  .< 

547.6 

510.4 

499.0 

462.4 

435.6 

409.6 

384.4 

360. C 

336.4 

313.6 

291.6 

270.4 

25Q.0 

230.4 

211  ^6 

193.6 

176.4 

160.0 

144.4 

129.6 

115.6 

102.4 

90.0 

78.4 

67.6 

57.6 

48.4 

40.0 

32.4 

25.6 

19.6 

14.4 

10.0 

6.4 

3.6 

1.6 

0.4 

O.C 

0.4 

1.6 

3.6 

6.4 

10.0 

14.4 

19.6 

25.6 

32.4 

40.0 

40.4 

57.6 

67.6 

79.4 

90.0 

102.4 

115.6 

129.6 

144.4 

160. C 

176.4 

193.6 

211.6 

233.4 

250.) 

270.4 

291.6 

313.6 

336.4 

360.0 

304.4 

4C9.6 

435.6 

462.4 

490.0 

518.4 

547.6 

577.6 

600.4 

640.0 

672.4 

705.6 

739.6 

774.4 

810.0 

846.4 

383.6 

921.6 

960.4 

1000.0 

STOP 

START 

»CNH«FT99F0l>l 

0.0 

99.0 

196.0 

291.0 

384.0 

475.0 

564.0 

651.0 

736.0 

019.0 

900.0 

979.0 

1056.0 

1131.0 

1204.0 

1275.0 

1344.0 

1411.0 

1476.0 

1539.0 

1600.0 

1659.0 

1716.0 

1771.0 

1824.0 

1875.0 

1924.0 

1971.0 

2016.0 

2059.0 

2100.0 

2139.0 

2176.0 

2211.0 

2244.0 

2275.0 

2304.0 

2331.0 

2356.0 

2379.0 

2400.0 

2419.0 

2436.) 

2451.3 

2464.0 

2475.0 

2404.0 

2491.0 

2496.0 

2499.0 

2500.0 

2499.0 

2496.9 

2491.0 

2484.0 

2475.0 

2464.0 

2451.0 

2436.0 

2419.0 

2400. C 

2379.0 

2356.0 

2331.0 

2304.0 

2275.0 

2244.0 

2211.0 

2176.0 

2139.0 

2100.0 

2059.0 

2016.) 

1971.) 

1924.) 

1875.0 

1824.0 

1771.0 

1716.0 

1659.0 

1600  .C 

1539.0 

1476.0 

1411.0 

1344.0 

1275.0 

1204.0 

1131.0 

1056.0 

979.0 

900.0 

819.0 

736.0 

651.0 

564.0 

475.0 

384.0 

291.0 

196.0 

99.0 

0.0 

STOP 
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M.  L.  BLODGETT 


HBUI  PM  CALCULATINS  SINGLE  MAT  FATA  STATISTICS  ASSURING  TNC  GAAIETT-IUNI  NOOEL  OF  INTERNAL  HATES 


IAFUT  ACTUS 


IE  Ft  AC  NCI  SOUNO  SFEEO  CCIEAO)  - 1SSA.S0  A/SEC 
AC RUSTIC  FRCaUEHCT  (FIE*)  - !••.•••  NEATS 

AAS  A AMCC  CHAD  ■ ISMS*  A 

NHAAEI  AF  AAHCC  STEFS  (AN)  - ISS 

TCEAA  OEFTN  (OEFTN)  - MAS.  A 

SUMACS  OEFTN  UKI)  - ISSS.S  A 

A ECS ITS A OEFTN  (IREC)  • ISSS.S  A 

IMF  NT  INTERNAL  HATE  FARAAETERS 


STRATIFICATION  SCALE  (0) 

SITAAFOLATSO  O-T  FAESUEHCT  (0 VIE 10) 
ESTRAFULATEO  FRAC.  FLUCTUATION  COELC) 
INERTIAL  FAEOUENCT  (FRCAIH) 

JSTAR  (JSTAA) 

CALCULATES  AUANTITISS 

RAASE  INCA.  (OS LA)  • ISO  .M  R 
OATENUROEA  CUT)  - S.AIOO  /A 

LIBAO  - ITS.SST  A 

ATS  OF  I/O  CATE)  • S.MSI 


RESULTS 

RAS  FNASS  FLUCTUATION  C FA I TOT) 

RAS  FNASB  RATS  FLUCTUATION  COOTTOT) 
RRS  INTENSITY  FLUCTUATION  (IOTOT) 

ERROR  FOR  FNASE  CSFNI)  ■ ( 

ERROR  FOR  FNASS  RATE  CEOOD  • I 

ERROR  FOR  INTSNSITT  (HOT)  ■ ) 


ISOS.  A 
S.SS  CFN 
SiME-S) 
S.SAES  CFN 

s 


S.ISSAS  CYCLE 
S. 141)2  CFN 
I.SOA  OR 


I NOS  I 

RAN6E-A 

MTCD-A 
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I(1SA)I(I)-A 

F HKD— CYCLE 
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1 

0. 

ISSS.S 

s.s 
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S.SSSM 

IS 

TOO. 

4T2.4 

-S.S2SSS 
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S.SS SSI 

S. SIS  12 

so 
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SS4.4 

•S.24ASA 

ISSS.S 

0.00*00 

S.SSSM 

so 

2 TOO* 

ITA.4 

•S.ITSSS 

SSST.S 

S.SIITT 

S.SS2SS 

M 

STOS. 

4S.4 

-S.STIT4 

2STT.S 

S.S20TS 

S. STAYS 

*0 

4 TOO. 

S.4 

-S.S1IS0 

24TT.S 

S.SAI2I 

S.IISTA 

40 

STOS. 

SI  .4 

S.SATTS 

24IT.S 

S.STTTS 

S.ISTSO 

TO 

ATM. 

144.4 

S.I44TS 

2IST.S 

S.STTSS 

S.1SSTS 

SO 

TMS. 

SSA.4 

S.I241T 

14ST.S 

S.lOSSS 

S.ISTSA 

TO 

STOS. 

MS.4 

S.2TSTS 

TTT.S 
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S.lOSSS 

too 

TMS. 

TM.4 

S.STSI2 
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S.100AS 

0.14120 

101 
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S.0 
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Fig.  C5— Plot  of  the  rmt  phaae  fluctuation  v*  horizontal  path  length 


PHIDDT(CHP) 


